Energy Dependence of Air Fluorescence Yield measured by AIRFLY 
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Abstract 

In the fluorescence detection of ultra high energy (> 10^* eV) cosmic rays, the number of emitted fluorescence photons 
is assumed to be proportional to the energy deposited in air by shower particles. We have performed measurements of 
the fluorescence yield in atmospheric gases excited by electrons over energies ranging from keV to hundreds of MeV in 
several accelerators. We found that within the measured energy ranges the proportionality holds at the level of few %. 
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1. Introduction 

The detection of ultra high energy (> lO^^eV) 
cosmic rays using nitrogen fluorescence emission 
induced by extensive air showers (EAS) is a wcU 
estabhshed technique [1]. Atmospheric nitrogen 
molecules, excited by EAS charged particles (mainly 
e^), emit fluorescence light in the w 300-400 nm 
range. The fluorescence detection of UHECR is 
based on the assumption that the number of fluo- 
rescence photons of wavelength A emitted at a given 
stage of a cosmic ray shower development, i.e. at a 
given altitude h in the atmosphere, is proportional 
to the energy 

^dep"""'(^) deposited by the shower 
particles in the air volume [2] : 

TVf "--(/i) = Et;^^-{h)Yair{X,Po,To)F{\,p,T), (1) 

where Yair {X,po,To) is the absolute yield (in number 
of photons per MeV) at a reference pressure po and 
temperature Tq, F{\,p,T) accounts for quenching 
effects, and p and T are the air pressure and temper- 
ature at the altitude h. Since a typical cosmic ray 
shower extends up to about 15 km altitude, the flu- 
orescence yield must be known over a wide range of 
air pressure and temperature. Measurements of the 
fluorescence yield dependence on atmospheric pa- 
rameters {F{X,p,T)) by AIRFLY are presented in 
a separate contribution [2,3,4]. 

Simple considerations suggest that fluorescence 
emission should indeed be proportional to the en- 
ergy deposited. In fact, the cross sections for elec- 
tron excitation of the 2P and IN nitrogen systems, 
which are the most relevant in the 300-400 nm range, 
are peaked at very low energies (tens of eV) and de- 
crease rapidly with energy of the electron (« E~'^ 
for the 2P and « logE/E for the IN). Therefore the 
fluorescence light induced by a high energy electron 
(> keV) will be mainly produced by the secondary 
electrons of eV energies. Since the total number of 
secondary electrons produced by the passage of the 
primary electron in the air volume is roughly propor- 
tional to the energy deposited, the fluorescence light 
is also expected to be proportional to the energy de- 
posited. The constant of proportionality should not 
depend on the primary electron energy. 

In fact, E^Jl"'^"'' {h) in Eq. (1) is the sum of the en- 
ergies deposited by EAS particles with a spectrum 
spanning from keV to GeV. It is thus important to 
verify the proportionality of the fluorescence emis- 
sion to the energy deposit over a wide range of elec- 
tron energies. Available measurements are limited 



to a few energies [6] or used indirect methods [7]. 
The AIRFLY (AIR Fluorescence Yield) collabora- 
tion has performed measurements of the energy de- 
pendence of the fluorescence yield at several accel- 
erators covering a range of electron kinetic energy 
from keV to hundreds of MeV. Results of these stud- 
ies are reported in the following. 

2. Electron energies from 3 to 15 MeV 

Measurements in the energy range from 3 to 15 
MeV were performed at the Argonne Wakefield Ac- 
celerator (AWA), located at the Argonne National 
Laboratory. The LINAC was operated at 5 Hz, with 
bunc;lies of maximum charge of 1 nC and length 
15 ps (FWHM) and typical energy spread of ± 0.3 
MeV at 14 MeV. The electrons exited the accelera- 
tor vacuum through a 0.13 mm thick beryllium win- 
dow. The beam spot size was typically 5 mm di- 
ameter, with negligible beam motion. The beam in- 
tensity was monitored with an integrating current 
transformer (ICT), immediately before the beam 
exit flange. The signal from the ICT was integrated, 
digitized, and recorded for each beam bunch. Fluo- 
rescence light produced by excitation of ambient air 
outside the beam exit was detected by a photomul- 
tiplier tube (Hamamatsu H7195 model) with a nar- 
row band 337 nm filter, located about 80 cm away 
from the beam axis. A shutter installed in front of 
the PMT allowed measurements of background. The 
PMT was surrounded by considerable lead shielding 
to reduce beam-related backgrounds. The accelera- 
tor timing signal was used to produce the integrat- 
ing gate of 200 ns width. Signals were recorded using 
a VME standard data acquisition system. 

The LINAC was operated in a mode allowing the 
bunch charge to fluctuate over a wide range. The cor- 
relation of the PMT and ICT signals, which showed 
a linear relation, was fitted and the slope Smeas was 
taken as an estimator of the fluorescence signal. The 
same procedure was applied with the shutter closed 
to estimate the background, which was subtracted . 

The measured fluorescence signal Smeas as a func- 
tion of kinetic energy is shown in Fig. 1. The full 
line is the expected fluorescence signal, Ssirn, es- 
timated by performing a full GEANT4 simulation 
of the experiment. In the simulation, the fluores- 
cence emission was taken to be proportional to the 
energy deposited by the particles in the gas. No- 
tice that the relativistic rise of the ionization losses 
in this energy range can be clearly seen thanks to 
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Fig. 1. Fluorescence signal as a function of kinetic energy, 
the full line is the result of a GEANT4 simulation where the 
fluorescence emission was proportional to the energy deposit. 
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Fig. 2. Relative difference between the measured and simu- 
lated fluorescence signal as a function of kinetic energy; open 
dots VdG data, closed dots AWA data. 

the accuracy of our data. The relative difference be- 
tween the measured and simulated fluorescence sig- 
nal, {Smeas " Ssim)/Ssim, is shoWU aS a fuUCtioU of 

energy in Fig. 2. The agreement between data and 
the Monte Carlo simulation confirms the propor- 
tionality of the fluorescence emission to the energy 
deposit between 3 and 15 MeV to a level of few %. 

3. Electron energies from 0.5 to 3 MeV 

Measurements were extended down to the mini- 
mum ionizing energy range at the Chemistry Divi- 
sion electron Van de Graaff (VdG) accelerator, also 
at the Argonne National Laboratory. The Van de 
Graaff accelerator was operated in pulsed mode at 
60 Hz, with beam currents from 0.2 to 0.8 /xA, and 
nominal beam kinetic energy ranging from 0.5 MeV 



to 3.0 MeV. The electrons exited the accelerator 
vacuum through a 0.152 mm thick dura-aluminum 
window. The beam spot size was typically 6 mm di- 
ameter, and a side-to-sidc beam motion of approx- 
imately 5 mm was observed due to small (< 1%) 
variations in the VdG energy on time scales of sec- 
onds. Fluorescence light produced by excitation of 
ambient air outside the beam exit was detected by a 
PMT located about 60 cm away from the beam axis. 
The PMT, shutter, 337 nm filter and data acquisi- 
tion system were the same as in the AWA LINAC. 
The beam intensity was monitored with the ICT de- 
scribed before and a Faraday cup. The total charge 
in the PMT was taken as a estimator of the fluo- 
rescence signal. To remove beam fluctuations, the 
PMT charge was normalized using the ICT signal. 
Background runs were also taken and substracted 
to the signal. 

A full GEANT4 simulation of the experiment with 
the Van de Graaf set-up was performed, and for 
each energy the predicted fluorescence signal Sgim 
assuming proportionality to the energy deposit was 
calculated. The relative difference between the mea- 
sured and simulated fluorescence signal, {Smeas — 
S Sim) /S Sim, is shown as a function of energy in Fig. 
2, together with the measurements of the AWA facil- 
ity. Notice that since measurements were performed 
at 3 MeV in both facilities, data are consistent with 
the proportionality of the fluorescence yield to the 
energy deposit with the same proportionality con- 
stant in the range 0.5 to 15 MeV. 



4. Electron energies from 50 to 420 MeV 

Measurements in the energy region of hundreds 
MeV were performed at the BFT (Beam Test Fa- 
cility) of the INFN Laboratori Nazionali di Fras- 
cati, which can deliver 50 to 800 MeV electrons and 
50 to 550 MeV positrons with intensity from single 
particle up to 10^ particles per bunch at a repeti- 
tion rate of 50 Hz. The typical pulse duration was 
10 ns. The beam exited the vacuum pipe through 
a 0.5 mm beryllium window, and produced fluores- 
cence light inside an aluminum pressure chamber 
(for a detailed description of the chamber see [2]). 
Given the low intensity of the beam (a few 10"^ elec- 
trons/bunch), a hybrid photodiode (HPD) with very 
good single photoelectron resolution was used to de- 
tect the fluorescence light. A 337 nm interference 
filter was placed in front of the HPD, together with 
a shutter that could stop the light for background 
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Fig. 3. Relative difference between the measured and simu- 
lated fluorescence signal as a function of kinetic energy. 

measurements. The beam intensity was monitored 
by Nal(Tl) calorimeter with excellent single electron 
resolution, placed at the end of the beam line. A 
fast scintillator was also used to monitor the beam 
intensity. The d(^pcndcncc of fluorescence light on 
the primary particle energy was measured in pure 
nitrogen in the range 50 to 420 MeV. We used ni- 
trogen to increase the fluorescence light yield, given 
the low beam intensity. The beam multiplicity was 
kept approximately constant at the individual en- 
ergy points. The fluorescence signal Smeas was esti- 
mated from the number of photoelectrons measured 
with the HPD, after background subtraction and 
normalization for the beam intensity. The relative 
difference between the measured and simulated flu- 
orescence signal, {Smeas " Ssim,)/Ssim, IS showu as 
a function of energy in Fig. 3, where Sgim is the ex- 
pected signal estimated by a GEANT4 simulation of 
the BTF set-up with the assumption of proportion- 
ality to the energy deposit. The agreement between 
data and the Monte Carlo simulation confirms the 
proportionality of the fluorescence emission to the 
energy deposit between 50 and 420 MeV to a level 
of few %. 
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Fig. 4. Relative difference between the measured and simu- 
lated fluorescence signal as a function of X-ray energy. 

trons of energies between 6 to 30 keV produced 
with this method deposit all their energy in a few 
mm of air. The fluorescence light induced by these 
electrons in the ambient air was detected by the 
photomultiplier, 337 nm filter and shutter system 
previously described, placed at 9 cm distance from 
the beam axis. The average charge recorded by 
PMT was taken as an estimator of the fiuorescence 
signal, after background subtraction. The X-ray 
beam intensity was monitored by ionization cham- 
bers placed along the beam axis. A full GEANT4 
simulation of the set-up, including the ionization 
chambers, was performed. The relative difference 
between the measured and simulated fiuorescence 
signal, {Smeas - S Sim) /S Sim, IS shown as a function 
of the X-ray energy in Fig. 4. Both for data and 
simulation, the fiuorescence signal was normalized 
to the ionization chamber signal. There is very good 
agreement between data and simulation, assessing 
the proportionality of the fluorescence emission to 
the energy deposit between 6 and 30 keV to a level 
of few %. 

6. Conclusions 



5. X-rays from 6 to 30 keV 

Fluorescence measurements with keV electrons 
were performed at the Advanced Photon Source 
(APS) of the Argonne National Laboratory. The 
intense synchrotron X-ray beam of the APS 15-ID 
line, after exiting the vacuum beam pipe to enter the 
experimental hall, produced an almost monochro- 
matic beam of electrons through photoelectric and 
Compton interactions with the ambient air. Elec- 



We presented measurements of the energy depen- 
dence of the fluorescence yield performed at several 
accelerators. We tested the proportionality of the 
fluorescence light to the energy deposited at a level 
of few % over the energy ranges 0.5 to 15 MeV, 50 to 
420 MeV and 6 to 30 keV. Notice that we performed 
only relative measurements within each range, and 
absolute measurements of the fiuorescence yield are 
in principle needed to verify that the proportionality 
constant is the same in the three measured energy 
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ranges [4] . Work in this direction is ongoing. On the 
other hand, given that the basic mechanism for the 
fluorescence yield is excitation by very low energy 
secondary electrons, it is hard to find any physical 
mechanism which could change the proportionality 
constant between 15 and 50 MeV. The AIRFLY data 
presented here would then indicate that the fluores- 
cence yield is indeed proportional to the energy de- 
posit for electron energies at least between 0.5 and 
420 MeV. Most of the EAS energy is deposited by 
shower particles within this energy range. 
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